Instabilities of monolayered epithelia: shape and structure of villi and crypts 
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We study theoretically the shapes of a dividing epithelial monolayer of cells lying on top of an 
elastic stroma. The negative tension created by cell division provokes a buckling instability at a 
finite wave vector leading to the formation of periodic arrays of villi and crypts. The instability 
is similar to the buckling of a metallic plate under compression. We use the results to rationalize 
the various structures of the intestinal lining observed in vivo. Taking into account the coupling 
between cell division and local curvature, we obtain different patterns of villi and crypts, which 
could explain the different morphologies of the small intestine and the colon. 

PACS numbers: xx 



An essential property of living tissues is a permanent 
cell turnover due to division and apoptosis, which has 
important effects on their mechanical response. Since 
cells often grow and divide in a constrained environment, 
cell division and cell death induce internal stresses in tis- 
sues that influence deeply their architectures and mor- 
phologies. Conversely, there is a mechanical feedback on 
cellular growth, differentiation, and organ development, 
which is nowadays an active area of research [1-6]. 

The intestine is the body's fastest renewing organ, and 
it is therefore a particularly interesting example to con- 
sider. It exhibits a variety of folded multicellular struc- 
tures [ ], called villi, which play a crucial role in favoring 
the exchange of nutrients. Although the development of 
the intestine has been studied for various animals, the 
formation of villi and their structure has not yet been 
described quantitatively [ ]. From a physicist point of 
view, these folded shapes in a periodic arrangement are 
strongly reminiscent of the patterns observed after the 
buckling of metallic plates. Euler buckling is the insta- 
bility leading to the lateral deflection of an elastic beam 
or a surface under load. A similar buckling instability 
is expected in growing constrained systems [ ], and has 
been invoked to explain fingerprints [ ] or the shapes of 
algae [ ] . Of particular interest is the recent report [ ] 
that intestinal structures can be reproduced in vitro with 
minimal physiological environment, reinforcing the idea 
that mechanical forces play a crucial role. Still, studies 
of cell renewal in the intestine take villi shapes as a given 
phenomenological curve [ ]. We argue that it is in fact 
a consequence of cell renewal itself. 

We present here a theoretical model for the intestinal 
structure based on the buckling instability of monolay- 
ered epithelial cells. We show that a few minimal ingredi- 
ents are sufficient to capture many key features of the in- 
testinal architecture. This model reproduces the correct 
orders of magnitude, the various patterns of folded struc- 
tures observed in vivo and gives insight into the physio- 
logical distinction between small intestine and colon. 



The intestine is a relatively simple organ: it is covered 
by a single layer of epithelial cells, which lay on top of 
the thin and relatively stiff basement membrane. The 
tube is surrounded by a soft connective tissue called the 
stroma. We adopt here a three layer model for the intes- 
tine, sketched on Fig. 1. 
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FIG. 1. Schematic description of an epithelial monolayer on 
top of a basement membrane and an elastic stroma. Crude 
estimates for the values of the parameters are [14] [ ], h c = 
10/xra, h b = lfim, E s = 400Pa, E c = 10 4 Pa, E b = 10 3 Pa. 



The underlying stroma is considered as an elastic 
medium of height H and Young's modulus E s . In prac- 
tice, the thickness of this connective tissue can be consid- 
ered as infinite as it is larger than the typical wavelength 
of the villi A ~ 100 fim. The basement membrane is a thin 
elastic sheet of height hb ~ 1/im, with a bending modu- 
lus Kb and an elastic modulus Ej,. The cell monolayer, 
has a thickness h c « 10/xra, and an elastic modulus E c . 
Considering the cell density as a constant, its bending 
modulus is K c = \E c h? c , which is much larger than the 
bending modulus of the basement membrane. The elas- 
ticity of the cell monolayer is dominated by curvature, 
and has no stretching energy, since the cells respond to 
strain by dividing or undergoing apoptosis over long time 
scales. We also consider that the cell monolayer can glide 
on the basement membrane. 

We first give a simple argument based on energy con- 
siderations which gives good insight into the mechanism 
of the instability and allows us to estimate orders of mag- 
nitude. The buckling instability is driven by the cell 
layer. When cells divide or die in a tissue, they exert 



stresses on their environment. The characteristic pres- 
sure is the homeostatic pressure Ph [ ], which is exerted 
on its surroundings by the tissue in its steady state, i.e., 
when cell division balances cell apoptosis. Since the in- 
testinal epithelium is a monolayer, the forces in the plane 
of the monolayer can be described in terms of a negative 
surface tension —7 related to the homeostatic pressure 
by 7 = Phh c . At the lowest order, the elastic energy of 
an undulation w(x) = wocos(qx) of the cell monolayer 
reads : 
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Buckling occurs when the energy of a finite wave vec- 
tor undulation becomes negative i.e. if the tension 7 
is larger than the critical value 7 C = (SK C E^) 3. The 
wavelength of the unstable undulation at threshold is 

then A c = 27r ( ^ J « 90/im, close to the observed 

in vivo wavelength (cf. 1). Remarkably, analytical so- 
lutions obtained by minimizing the full elastic energy 
above threshold, for an arbitrary pressure, retain the 
same wavelength. The critical homeostatic pressure for 
villi formation is typically Ph c = 1400 Pa.m. This value 
is compatible with a homeostatic pressure in the range 
10 3 — 10 4 Pa [16]. One can therefore expect the intestinal 
lining to buckle under the pressure of the epithelium. 

Surprisingly, animals of very different sizes (mice and 
humans for instance), have villi of equivalent dimensions. 
This makes sense in our description, the size and wave- 
length are only dictated by elastic considerations, and the 
elastic parameters are roughly constant in most mam- 
mals. 

The previous energetic approach reproduces accurately 
the undulations of the small intestine. The architecture 
of the colon is visually very different [ ]: there are no 
villi but only crypts extending into the stroma. 

This asymmetry prompted us to take into account the 
non-uniform division rates along the villi. In vivo, cells 
multiply in crypts from stem cells and undergo apoptosis 
at the tip of the villi. They constantly flow from the 
crypts to the villi, differentiating along the way. The time 
scale for complete renewal of an epithelium is 3 days in 
the intestine, and 7 days in the colon. We do not wish to 
enter here in the molecular details, we will suppose that 
the non-uniform cell division can by modeled by a growth 
rate depending on the curvature of the cell monolayer, a 
hypothesis which has biological grounds [13]. 

We use as a reference state the homeostatic state 7 = 
70 where the monolayer is flat and reaches a steady state 
so that the cell division rate kd balances the cell apoptosis 
rate k a , then we expand our equations around this state. 
Since the monolayer is polar, the apical and basal sides 
are not equivalent and a linear coupling with curvature is 
allowed : k^ — k a = — £(7 — 70) + aV 2 w. The coefficients 
£ and a measure the dependence of the growth rate on 
pressure and on curvature respectively. They are positive 



so that the growth rate decreases with pressure and is 
larger in the crypts. Assuming that the cell density in the 
monolayer is essentially constant, the cell conservation 
equation reads: 



V.v = kd - k a 



-£(7 -7o) + a\/ 2 w 



(2) 



The pressure gradient in the monolayer is balanced by the 
friction force with the basement membrane: (v = — V7, 
where ( is a friction constant. These two equations 
give the pressure in the cell monolayer as a function 
of the undulation amplitude. For small friction coeffi- 
cients 7 = 70 — a(w while for large friction coefficients 
7 = 7o + f V 2 w. It is difficult to obtain precise values of 
the parameters £, £ and a and therefore to decide which 
regime is more relevant. We have performed calculations 
in both regimes, and observed the same physical behav- 
ior. In the following, we present calculations in the small 
friction regime, and define p = a(. 

For large values of a, the crypts and villi become asym- 
metric. Villi are flattened, since the pressure is locally 
lowered by apoptosis, and crypts are deepened, since the 
pressure is locally increased by division (Fig. 2 ). 
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FIG. 2. (a) Small intestine morphology, showing developed 
villi, and (b) colon morphology, showing only crypts. Units 
are micrometers. 

We start the analysis with a = 0, in order to study 
the morphology of the small intestine. A large a is then 
included to explain specificities of the colon. Our main 
result is a phase-diagram which presents the possible in 
vivo morphologies of the intestinal tube. 

For each point in the monolayer, we define a displace- 
ment vector u = (iii, ^2, w), where w is the vertical dis- 
placement of the interface and (111,11,2) the lateral dis- 
placement. We obtain the steady state shape and struc- 



ture of the villi from the force balance on the basement 
membrane. We introduce in the force balance a local 
drag force proportional to the local velocity of the mem- 
brane, which controls the relaxation of the membrane to- 
ward its equilibrium shape. The actual dissipative force 
is certainly non local and more complex but as we are 
interested only in the steady state shape, this provides 
a convenient way to relax, in our numerical solutions, 
the system toward its steady state configuration. The 
basement membrane is described by the classical Foppl- 
Von-Karman [19] equation, which applies to thin rigid 
plates at moderate deflections. We consider for simplic- 
ity a membrane with vanishing Poisson modulus. 

After rescaling all lengths by the thickness of the base- 
ment membrane hb and time, the equations of motion of 
the basement membrane are : 
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The various terms in equation (3) describe the fric- 
tion force, the curvature of the cell layer, the pressure 
exerted by the cells, the stretching force of the basement 
membrane, and the elastic force due to the stroma. The 
elastic force of the stroma £ e i is proportional to the mem- 
brane displacement u. Although we are not writing it for 
the sake of simplicity, we also took into account the non- 
linear dependance of the curvature and tensions forces 
on Vu>, deriving them from their full Hamiltonian. This 
does not change the qualitative results, but raises the 
issue of which non-linear elasticity should be chosen for 
the stroma. This is beyond the scope of this paper and 
should be addressed in future work. The linear relation 
for £ e i is best written in Fourier space for a wave vector k 
as f e / = 3§^(^)^Mu where M is the 3 x 3 matrix given 
by 

/ (2k 2 - k 2 y )/k k x k y /k \ 

k x k y /k (2k 2 -k 2 x )/k (4) 

\ 2k J 

Whereas, the buckled state is a simple sine-like func- 
tion in ID, there are a variety of possible patterns in 
two dimensions [20-22]. The three main morpholo- 
gies are finger-shaped villi, herringbone, or disorganized 
labyrinth patterns. Following [ ], we implemented a 
semi-implicit integration method. All products are done 
in real space, and all derivations in Fourier space, for 
maximal efficiency. As shown in Fig. 3, we obtain a 
phase diagram regrouping different equilibrium solutions 
for characteristic physiological values of our parameters. 
The parameter that we vary here is the buckling pressure 
exerted by the monolayer, on the horizontal axis. 

At low buckling pressure, finger shapes are stable, 
but for large stresses, it becomes difficult to bend the 




FIG. 3. Phase diagram showing the various possible villi 
organizations. Fig. A is colon-like whereas B, C and D are 
small-intestine-like, respectively finger-shaped, herringbones 
and labyrinth. The vertical axis is the coupling between cell 
division and curvature and the horizontal axis is the pressure 
exerted by the cell monolayer. 



membrane in all directions, and labyrinth or herringbone 
structures are favored. This has been observed experi- 
mentally for various inert thin film systems [ ]. And 
indeed, villi are not always finger-shaped. They some- 
times adopt complex folded patterns strikingly identical 
to the ones that we predict [ ] [17]. 

To quantify the buckling pressure, a good proxy is the 
height of villi, as indicated by equation (6). Although 
villi wavelength is roughly constant throughout the intes- 
tine, villi height is not. At the entrance of the intestine 
(near the stomach), villi are very long, but their size de- 
creases gradually in more distal parts [ ] . This leads us 
to conclude pressure decreases along the intestine, justi- 
fying the relevance of our phase diagram. 

As an example, Fig. 4 shows a transition between 
finger-shaped and herringbone villi, in vivo [17], and from 
our numerical integration when increasing the pressure 
along the horizontal axis. 
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FIG. 4. Transition between finger and herringbone villi, com- 
paring a transverse cut (left) and our simulation (right). Pres- 
sure increases along the horizontal axis from left to right. 

We therefore expect that labyrinth-like villi are more 
likely to be found near the stomach, and that in distal 
parts, villi should be finger-shaped. Indeed, this is what 
physiological studies show [ ]. 

While we do not claim that this is the mechanism for 
the formation of villi during embryogenesis, we argue that 
when the system reaches a steady state, stresses must be 
balanced, as described by our equations. In fact, several 
experiments have shown that when increasing cell apop- 
tosis for a short time [25], villi disappear, and in a matter 



of days, regain their initial equilibrium state. 

Seeing villi as a product of a buckling instability could 
shed new light on intestinal pathologies. It should be 
noted that in our model, the elasticity of the basement 
membrane plays a crucial role, since it stabilizes the am- 
plitude of the folds. In coeliac diseases, atrophy of the 
villi correlates with an increase of the apoptotic index, 
but also with an increase of the basement membrane 
thickness [ ]. This translates into a lower pressure, and 
a higher energetic cost to bend the membrane, which 
can cause the pressure to decrease below the buckling 
threshold, and thus the abrupt disappearance of villi. 
Conversely, if the basement membrane breaks, which is a 
necessary event in colon cancer, this model predicts un- 
controlled outgrowth of the intestinal lining, similar to 
polyp formation 

In addition to this equilibrium state, there are sec- 
ondary buckling events in crypts that result in a much 
higher crypt to villi size ratio than described here. This 
mechanism has been investigated in Ref.[26] from a phys- 
ical point of view and we will ignore it in this work. 

The colon can be studied by considering a high asym- 
metry in cell division rate due to curvature. The pres- 
sure is then a function of the deformation w as 7 ~ 
70 — pw(x, y, t). The epithelium renewal rate in the colon 
is about twice slower than in the small intestine, sug- 
gesting that the stresses exerted are also lower. Also 
relevant is the fact that CDX1 and CDX2 levels are 
higher in the colon [ 1 8] , and that these protein inhibit 
j3 - catenin, thus favoring earlier differentiation. Physi- 
cally, this would mean that high stress regions are more 
localized in colonic crypts than intestinal crypts, and that 
the asymmetry parameter p is higher in the colon. 

Integration of these new equations for several values 
of p demonstrate, as shown Fig. 2 and Fig. 3 A and 
B, a transition between intestinal and colonic morpholo- 
gies. This transition is sharp, as observed in vivo at the 
intestinal-colonic junction. This is consistent with the 
following analytical calculation made for small values of 
pin a one-dimensional geometry. 

In one dimension the monolayer deformation ampli- 
tude follows the equation 

= K c A 2 z + joAz - E b h b V.{{Vz) 2 Vz) + f el - pz(Az) (5) 

For small values of p, we expand the solution as a sum of 
two cosine functions: h(x) = ho + hicos(qx) + h2cos(2qx), 
with q — ^ and identify the terms in cos(qx) and 
cos(2qx). The amplitudes are such that h\ <C h\ and 
pli2hi <C h\. Within this approximation, the value of h\ 
is unchanged compared to the symmetric p = case : 



larger cos(x) term, the troughs are deepened and the 
crests flattened. 
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Our numerical solutions predict that intestinal villi are 
organized in a square lattice (as can be proven analyt- 
ically by energy considerations, following [22]), whereas 
colonic crypts are organized on an hexagonal lattice as 
seen on Fig. 3 (in agreement with weak crystalliza- 
tion theories [ ]). Although the system is noisy, this is 
strongly hinted in vivo by the fact that transverse cuts of 
the small intestine do show square-shaped villi, whereas 
crypts are round and have a higher number of neighbors. 

This simple buckling model predicts a full phase dia- 
gram explaining most of the structures observed in the 
small and large intestines. This suggests a deep con- 
nection between tissue architecture and the stresses pro- 
duced by the dividing cells. Intestinal shape and renewal 
are at least partially controlled by a mechanical balance, 
which is disrupted in the case of intestinal diseases. Since 
it has been proven that an excess of mechanical pressure 
could, in itself, induce colonic cancer [ ], a detailed un- 
derstanding of this balance is of great interest for cancer 
research. 

We thank S. Fre and M. Huygue for discussion and for 
showing us in vivo samples, and A. Berges for help with 
the manuscript. 



The non-symmetric term fi2 is negative, in agreement 
with what is expected: when subtracting cos(2x) to a 
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